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ABSTRACT 


To assess criteria for the determination of actual evaporation 
by remotely sensed surface temperature, a sensitivity analysis of the 
energy balance equation is performed. 

Constraints in the design of field experiments are specified and 
criteria to collect direct spot measurements of actual evaporation 
are formulated. 

By Taylor expansion of the energy balance equation, analytical 
approximations are derived and fitted to experimental data. The latter 
are obtained by the Bowen ratio for both actual bare soil evaporation 
and actual transpiration from natural vegetation. 

A comparison of different equations, with actual evaporation as 
a function only of surface albedo and surface temperature is also 
presented. 


I. INTRODUCTION 


The work contained in this note was performed while the author 
was on leave from TECNECO spA (Italy) for eight months, during 1979. 

According to many recent references, remotely sensed surface 
temperature allows for a rather promising way. to the estimation of 
actual evaporation losses from large non-homogeneous surfaces. However, 
plots of direct measurements of actual evaporation versus surface 
temperature show a substantial scatter of the points. Therefore it is 
not easy to define some analytical expression of the relationship be- 
tween evaporation rate and surface temperature. The problem is that 
such a scatter is due to the interference of other independent para- 
meters of the surface energy balance and due to the lack of some 
a-priori equation to relate the parameters involved. 

A way to solve this problem is to derive analytical expressions 
by Taylor’ s expansion of the energy balance equation. 

II. TAYLOR'S EXPANSION OF THE SURFACE ENERGY BALANCE EQUATION 
a. Concept 

According to the previous suggestions, it seems worthwhile to 
specify in a formal the concept of relationships between evapora- 
tion rate and the different variables of the surface energy balance 
equation. There is however no need to establish a new relationship 
as it is already implicitly included in the energy balance equation 
at ground level. 


I 



p Q 

LE • CT a -T g ) ♦ ( l-(i)K g ♦ t:*OT a - eoT g •* G (.1. cm" 2 . day” ' ) (1) 

where L * latent heat of vaporization of water (at mean air tempe- 
rature) (J.g” 1 ) 

E • evaporation rate (g.cm .day ) 

P a c g ■ thermal capacity of air (J.cm” 3 .°C~ 1 ) 
r g « surface aerodynamic resistance (day. cm” 1 ) 

- air temperature (K) 

T g ■ surface temperature (K) 

a ■ surface reflectivity 

R g * shortwave incoming radiation (J.cm .day ) 
e* * air apparent emissivity 
e ■ surface emissivity 

o ■ Stephan- Boltzmann constant (J.cm”^,day -1 . 0 K”S 

• - 2-1 
G ■ heat soil flux (J.cm .day ) 

Such an equation may be thought as a hypersurface in with 
the coordinate axes being E, T , r , T_ , a, R , G. Eq. (1) defines 

(tab S 

a general relationship between evaporation rate and the meteorological 
conditions at the surface level. This relationship is not subjected to 
any kind of constraints related to the kind of surface and/or meteo- 
rological situation. That is to say if the turbulent exchange coeffi- 
cients are evaluated without restricting hypothesis. The obvious dis- 
advantage is the extensive data requirement of this equation: compre- 
hensive and area-spread meteorological data are needed for the evalu- 
ation of regional evaporation losses by eq. (1). To evaluate regional 
evaporation, using, if possible, only remotely sensed data over large 
areas we thus have, for practical reasons, to look for more simple re- 
lations. This is equivalent with defining functions of a degree of 

dimensions lower than eq. (1), for instance a line-relationship (E, T ), 

s 

or a surface-relationship like (E, T g , rt). These functions can be 

defined through a Taylor's expansion of eq. (1) around some physical 

status E*, where E* r f(T*, a*, r*, T*, C * , R* ) with a notation in 
1*7 s a a s 

• 

Taylor's theorem, using the vector £ as representing the point 
with coordinates (x,, inIR cai be written as: 
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f(s) - f<$*) ♦ f * (s*) (s^B*) ♦ f"<S*)<S-«.*) /2 J ♦ f*V" ^ (s*) (8**8*)^"*^ / 

(n- I ) l * R n (2) 

where: R is the remainder 
& 

t 1 is the total i-th order derivative with respect to s 5 x^» 

..... \ 

The sum of an unlimited number of terms as in eq. (2) approximates 
the function f if, 'and only if, the remainder approaches zero as the 
number of terms becomes infinite. Such a series is convergent to the' 
function f with the number of terms. Then the error in the approxima- 
tion of the function f depends on such a number. 

A Ttylor's series can be used to approximate the function LE * 

f(T , a, r i Ta, G, R ). For oui purposes it is very convenient to 
s a s 

use only the linear terms in eq. (2), i.e. applying only the first 
order derivatives. Then from eq. (2) we get: 

E » E* + D E dT i ,i + D E da i,.i 7 + D E d t, i_.i- + 

+ °Ta E dT a V±7 + °G E d G — 5*—7 + °Rs E d R s V±7 (3) 

where: E* is the value of eq. (1) at (T *, a*, r *, T *, G*, R *) 

s a a 8 

D . is the partial derivative with respect to x. 
xi i 

i. is the versor (unit vector) of axis x. 

-J J 

i-.i^, is the scalar product between the two versors (then 

3 • A 

equal to the cos. of the angle i. i.) 

“2 “1 J * 

E, E* are in g.cm .day ; i_ 7 is the versor of E 

b. Numerical evaluation of “'the derivative 
with respect to the surface tempera- 
ture 

The geometrical framework info^ can be followed further to obtain 
a formal expression of the relationship between E and T g : 

-2 - J 

E(T g ) - E.Ug.ij) (g.cm .day ) (4) 
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E(T') it the projection of the vector & on the coordinate plane 
i_, ij , Such e projection can be approxiaated according to aq. (1), 
(3)* 

E* (T ) - E'* - * 4 e 0 T «> dT . (g.cm -2 . day" 1 ) (5) 

3 la r • ■ 

s 

where: E‘(T ) is the first-order approximation of E(T ) and 

s s 

E'* the value with T g - 

This kind of equation can be applied to fit experimental data on 
evaporation rate vs. surface temperature* As can be seen from eq. (5) 
the slope of the function E'(T ) depends on T (supposed to be known 
everywhere) and the prevailing meteorological conditions over a wide 
range of the surface roughness. 

However, in the conaon range of surface temperatures * low values 

of the aerodynamic resistance may cause a shift from a small slope of 

E'(T ) to a steeper one. In Table 1 values are reported for different 
s 

aerodynamic resistances, as calculated from: 

P c 

- j- ♦ 4 e o T^) (g.cm 2 .day *.K *) (6) 

L 1 S 

a 

E'(T g ) is decreasing everywhere, since quantities inside brackets 
are positive, when 6T g is positive. Anyhow it follows from Table 1 
chat experimental points with low values of r' require a steeper slope 
of eq. (5) in order to be fitted. 


Table 1. Values of the slope of E'(T g ) - eq. (6) - for different values 
of the aerodynamic resistance r 

a 


r a 


-D„ E 
Ts 

(day. cm'’ * ) 

(B.m -1 ) 

(g.cm ^.day *.K *) 

1.16 

(1000) 

.0212 

. 1 16 

(100) 

.0249 

.01 16 

(10) 

.0617 
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C, Variability of evaporation rata at 
related to the other meteorological 
variables 

In Fig- I data of actual soil evaporation as obtained by the 
Iowan ratio method are plotted versus the corresponding surface tempe- 
ratures. 


actual evaporation £ 
(mm Soy*) 
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Fig. I. Evaporation rates as calculated according to eq. (5) with a 
constant slope (middle broken line). The range of scatter due 
to surface albedo as calculated by eq. (7) is shown by the 
short broken lines. Heasured evaporation rates E are plotted 
versus measured averaged daily surface temperature T . Dotted 
points are within, open circles are outside the range bounded 
by eq. (7) 


The scatter due to differences in surface albedo is also shown, 
as evaluated from; 

R , 

D ft E - - ~ da (g, cm” 2 . day ) (7) 

and evaluated with 5a • (a x -a lJ )/2, corresponding to a 6 q E ■ ♦ 1.4 
wm.day” 1 , where a and a are the maximum (0.583) reap, the minimum 

X H 

value (0.262) of the surface reflectivity in the data set. 
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«clu*< «»a» c rot »on E 



lurfoc* t*mp*fotuf* T # 


Pig, 2* At Pig. 1, but the range of scatter is not only due to surface 
albedo but also due to air temperature as calculated by eq. 

( 8 ) 


In Fig. 2 the scatter due to the air temperature, according to 

the derivative t 


* * Z<* * T . 

A 


t ~ 2 „ -i\ 

(g.co .day ) 


is elso shown. In eq. (8) 6 T ■ (T x -T‘ N )/2, where T* and T N are the 

a a a a a 

maximum rasp, the minimum value of air temperature in the data set. 

As can be seen from Fig. 2, this may account for a variability in 
evaporation of an additional ♦ 0.3 mm. day A comparison of Fig. I 

with Pig* 2 shows that the scatter due to the variability of surface 

reflectivity could be justified by a curve E (T g ,a). 

Before continuing, the dependence of D Tg E on T & must be evaluated. 
In the Pigs. 3 and 4 E(T g ) is calculated and plotted with a varying 
slope (eq. 5), while the scatter is calculated according to: 

R X N »£' O f 3 I* - l" pc T* - T* 

( 4t) - *2-2-. _ 1 - r"“ (-W-1) 

MIX As A 


4 e' o T 3 


pc t - r 
x a a / a a x 

' " Lr ' 2 * 


(g.cm” 2 .day" 1 ) 


As can be seen, the results of Fig. 2 and Fig. 4 are practically 
equivalent. 
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Fig. 3. As Fig. I, but with the slope dependent on surface temperature 
T 

s 



Fig. 4. As Fig. 2, but with the slope dependent on surface temperature 


In a way similar to eq. (4), an analytical expression of E(T^,a) 
can be derived: 


E(T s ,ct) - E.(i 8 ,ij,i 2 ) 

or 


, " 2 j ~K 

(g.cm .day ) 


( 10 ) 
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(to 


E(T s ,a) - E* ♦ D Ti E dT # ♦ D ft E d a (g. cm" 2 , day' 1 ) 

d. Numerical calculation of th« deriva- 
tive with respect to surface tempera- 
ture and reflectivity 

Eq. (11) can be written as: 

I ^a £ a 3 _ o _ I 

E(T ,<x) - E* - £(-— d T s M E n f d I 4 ♦ R da) (g.cm -\day 1 (12) 

A 

This equation corresponds to eq. (6) and, according to the sane 
procedure applied there, the zeros of the second term of the right 
hand side of eq. (12) can be evaluated: 

Pc n 

-£-* dT ♦ 4 c fl T d T. ♦ R da-0 (J.co“ 2 .day” 1 ) (13) 

r s s s s 

;i 

Solutions T g " of eq. (13) depend on the combination of the incre- 
ments 6t and 5a. In Table 2, values T " are shown as included in the 
s s 

interval corresponding to £.he different values of r (Table 1). The 

fi 

left hand side of eq. (13) can be < 0 when either 5T g or 6a is negative, 
according to the value of r . The headings of rows and columns in 

a 

Table 2 show the different choices of 6T and 5a. 

8 

For a proper interpretation of these results it must be kept in 
mind that they have been derived by Teylor's expansion of eq. (I). 

Thus the trend of the evaporation rate £ is, strictly speaking, de- 
termined cily in the neighbourhood of the point one is dealing with. 

After collecting experimental values of the evaporation rate for 
different surfaces and periods one by one, the procedure to be applied 
is the following. Firstly one checks for a point Pj if ^Ej 4 
0 or >0. Secondly, for the points P. -*■ P, if 6_ e| <0 or >0. P l 

j ' ^ tSgOt p 

This means to check if the points P| and can be I connected by 

a raonotonic single-value function of the couples (I ,a). The points 

. 2 1 
p and F, constitute a data set denoted S . Then the same procedure 

1 ^ 2 
if applied to a certain point looking for the nearest point in S 

with respect to the coordinate system T g ,a. Now one has to check if 


8 


fable 2. Surface temperatures T f H (K) satisfying eq. (13) for different 
choices of $T g and fa, with 0.01 16 r g 4 1.16 (day.ca* 1 ). 

When 6T g < 0, E > E for T g > T g ". The reverse holds for 6T g > 0 



fa <• 

.1 

i 

a 

5 

.1 

K 

(°C) 

K 

(°C) 

6 T - 1 

s 

-553<T g M <-488 

(-826<T g "<-76l) 

396<T g w < 386 

( l23<T g "< 213) 

o> 

HI 

• 

B 

1 

396<T g * , < 486 

( 123<T g "< 213) 

-553<T g H <-488 

(-826<T g "<-761) 

6 T g - 5 

-424<T g "<-287 

(-697<T g ”< 560) 

-31 KT g ”< 283 

(-584^T g "< . 10) 

6 T s • -5 

-31 l<T g "< 283 

(-584<T g "< 10) 

-424<T g "<-287 

(-697<T g "< 560) 


2 

P. can be connected to its nearest point in S with a slope of the sene 

■* 2 2 3 

sign as in S . If so, add Pj to S to get S . If not, add P^ to a second 

data sub-set where the slope is of opposite sign. Then repeat the pro- 
le- ! 

cedure for each P^ with respect to S . Then it is possible to fit 
independently the two sets with positive (6 Tg a E>0) and negative 
(6 t> ft <0) variation. 

In the following, different relationships between E, T , a will 
be analytically derived and fitted to data obtained from a series of 
experiments performed in an inhomogeneous desert area. 


e. First-order approximation of the 
surface E(T g ,ot) 

In the present case the well known condi -.ion for four points be- 
longing to the same plane leads to the following expression of the 
linear function E'(T ,ct): 


1 

1 


V T . 

a-a 

E-E 

_1 2 
V T . 

1 2 

e’-e 

a -a 

T 3 -T? 

3 2 
a -a 

e 3 -e 


s s 


114 ) 
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The point. P, ? (T*, a 1 , B 1 ), Pg S (tJ, a 2 , B 2 ), P £ (tJ, a 3 , E 3 ) 
Must be chosen on the bssis of eq. (I). This sssns that, with a Minimal 
quantity of data, tha Figs, I to 4 togathar with eq. (I) can be used 
to choose tha points P ( , Pg, Pg. Tha terms B l^.i^ and E i^.ig in aq. 

(10) cm be represented on the coordinate planes (£, T § ) and (B, a) 
respectively. If s constant value for D^E is used, E'CT^) and E'(a) 
are straight lines in the two planes. Evaluation of the slope of E'(a) 
requires an estimation of the shortwave incoaing radiation. In this 
paper the average R # of the whole data sat (winter tine and Sumner 
time data) is used. Evaluation of the slope of E'(T s ) is note cumber- 
some, because * f As can be scan from Figs. 3 and 4, 

D E is not varying strongly in the interval 280<T <320 (K) (7<T <47 °C). 

Ts * s 

As far as r # is concerned, there is a strong variation of r # in the 

existing data set. This problem will be discussed later in detail in 

section TII.c. Anyhow, in order to evaluate the feasibility of very 

simple formulas (as derived in a formal way fre? •‘he energy balance 

equation) only mean values are used. 


III. NUMERICAL DEVELOPMENT OF LINEAR EVAPORATION EQUATIONS 

a. Selection of suitable experimental 
points 

Whether or not the use of eq. (14) is successful depends on the 

selection of the points P )f P^, (section II. e). The most reasonable 

choice seems to use extreme values either for T or a. One of the 

s 

points P. can be on the axis E, the other two on the coordinate planes 

(E, T # ) and (E, a). With this solution it is likely that the experimental 

values of E are closer to the plsne E'(T ,a). The first step is to 

s 

calculate the intercepts of the straight lines E* (T ) and E'(cx) from 

s 

eq. (S) and (7). When using the input data presented in Table 3, the 
following values for the derivatives are found: 

D_ E * -0.206 (tmn.day -1 .K~* ) D E * - H.8* (ma.dav’’ 1 ) (15) 

Ts a ‘ 
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Table 3. Numerical input data used for the 
E*(T a ) and E'(o) 

f- v • ■ • , 

calculation of the functions 

Parameter 

Value 

Units 

P c 

a a 

1.154 

-7 

J.cm" 3 .»f 1 
J.cm^.day" 1 .^ 4 

4 e a 

18.615.10 1 

*a 

2.15 

. -1 
day. cm 

T 

s 

300.8 

K 

V 

2154.8 

J.co .day 

L 

2432.3 

J.g” 1 (at 27°C) 


The functions B*(a)» E*(T g ) are described as: 

E»(T ) - -0.206 (T -280) + m (mm.day 1 ) (16a) 

E’(a) ■ -8.86 (cr.25) + n (mm.day *) ( 16b) 

The intercepts m and n can be calculated simply by imposing that 

tvo known points belong to the lines. The origin of the coordinate 

o 

planes was chosen corresponding to a * 0.25 and T # ** 280 K, both 

values being smaller than the minimum values in the data set. The point 

Qj S (280, 0.362, 3.8) is on the plane (E,a) and the point Q 2 = (300.8, 

.25, 3.8) is on the plane (E,T fi ). The average values are: I » 3.8, 

T * 300.8, s « 0.362, and accordingly m * 8. I and n ■ 5.0 were found. 

8 „ 

One of the three points required was chosen as E * (m+n)/2. Then an 

arbitrary choice of the points P(T , a, E) can be: 

s 

P J = (280., .25, 6.6); ? 2 = (300.8, .25, 4.); P 3 = (280, .7, 1.2) 

The point P 3 va * al *° calculated from eq. (16b) on the plane (E,a). 
From such a selection and according to eq. (14), a relationship between 

E, T and a was obtained as: 

s 

E'(T ,a) - - 12 . a - 0.125 T + 44.6 (mm.day l ) (17) 

8 8 

Since a poor agreement wai, obtained for experimental values cor- 


1 1 


responding Co high reflectivities and low temperatures, a new trial 

ve* 8on*withr ( - 

Pj 5 (280, .25, 6.6); ? 2 5 (300.8, .25, 4.); Pj = (282.9, .583, 3.4) 

where P ^ is measured. In this case the relationship obtained was: 

E'(T ,o) - -8.52a - 0.I25T ♦ 43.73 (mudey" 1 ) (18) 

8 8 

b. C o a p hr ison between calculated and 
experimental values 

The accuracy of eq. (18) and eq. (5) was tested against the expe- 
rimental points of the Figures 1 to 4. In Table 4 measured E-data, es- 
timated values B* (T ) from eq. (5) and E'(T .a) from eq. (18) are com- 
8 8 

pared. Standard deviations are also shown. The experimental values 
were collected over quite different surfaces and different periods of 
the year. 


Table 4. Comparison between experimental data with values obtained from 
eqs. (5) and (18) 


E-function 

Mean 

(mm. day *) 

0 

(mm. day *) 

E 

3.8 

1.26 

E' (T ) from eq. (5) 

4.1 

1.43 

E-E'(T) 

s 

-.3 

1.46 

E'(T g ,CX) from eq. (18) 

3.2 

1.21 

E a - E ’ (T s ,a) 

0.6 

1.1 


According to the theory developed in section II. d, it is still 

necessary to check the sign of the derivative for every point in the 

data set. Such a check was performed, with the derivatives evaluated 

using the average shortwave radiation K , and the actual data of T 

s s 

and r a » The aerodynamic resistance r g was evaluated according to 
FEDDES (1971) as: 


12 


r . ff(l) 1.15 u”* 7 ^J (s.m“ J ) (19) 

* p a L J 

• kImn e* * T«tio of 'molecular weight of water vapour to dry air 

p * atmospheric pressure (pascal) 

* “I 

u • wind Velocity (m.s ) 

f(l) - a l b 

1 « height of roughness elements 

a * 0. 167 x 10 7 when 1 * 20 cm 

a ■ 0.3704 x 10 7 when 1 > 20 cm 

b • .59 when 1 ^ 20 cm 

b ■ .2827 when 1 > 20 cm 

The increments 6T g and ht were evaluated according to the proce- 
dure described in section II. d. 

With such a procedure it was possible to separate the experimental 
points in two sub-sets, according to the negative and positive values 
of D_ E. The sub-set with the positive derivative included five 

18 *Cl 

points out of the 33 presented in Figs. 1 to 4. A new evaluation of 
the accuracy of eq. (18) was performed for the sub-set with the negative 
derivative. Results are shown in Table 5. 


Table 5. Comparison between experimental and calculated values of the 
sub-set with the negative derivatives 


E-function 

Mean 

(mm. day *) 

0 

(mm. day *) 

E 

3.6 

1.25 

E’(T g ,a) 

3.1 

1.2 

E g -E'(T g ,a) 

0.5 

l.l 


The results of Table 5 are not excellent, but they still support 
the theory. In the following sections the effect of air teng>erature 
and the way how to handle data with very low r , will be analysed. 
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Effectiveness 
sistance, r 

A 


o f 


the aerodynamic re- 


Up till now the dependence of E upon r. was not taken into account 
Such a hypothesis seems rather inaccurate because thus points of high 
aerodynamic resistance were connected with points of low aerodynamic 
resistance. From eq. (I) D r E is derived as: 


D E - -i PC (T -T ) 4 dr 

IT Li 8 A A S L d 

* r 

a 


(g. cm* -2 . day" 1 ) 


( 20 ) 


Since a finite evaluation of 6 E is needed, it is better to write 

r a 

eq. (20) as: 

p. Tp. 

J I J 


E 


J- p c (T -T ) 
L a a a s 


P. a 

l 


(g.cm 2 .day ') 


( 21 ) 


where P^, P^ are points in the data set, and the average (T a ~T g ) is 
taken between P^ and Pj. The correction according to eq. (21) must be 
evaluated for the points in the data sub-set where D_ < 0, for those 

is jtZ 

points showing the highest roughness. The points to be considered to- 
gether with the data required for the calculation are depicted in 

Table 6. 


Table 6. Experimental input data corresponding to days and sites with 
low aerodynamic resistances r . The last column represents 
data obtained from calculations with eq. (18) 



T 

s 

(K) 

T 

a 

(K) 

a 

r 

a 

(day. cm ) 

E 

(mm. day ) 

E'(T s ,a) 
(mm. day l ) 

1 

308.7 

305.7 

0.298 

0. 1 1 

5.5 

3.58 

2 

309.2 

312.7 

0.301 

0. 1 1 

6.3 

3.54 

3 

309.9 

310.3 

0.303 

0.13 

6.4 

3.50 

4 

307.9 

306.2 

0.304 

0.07 

5.2 

3.56 


14 


Until now the relationship between actual evaporation rate E and 

surface characteriatics T and a has been sought in a three-dimensional 

...... 8 *• * - ' * 

space (E , T , a). In principle evaluation of the effectiveness of r 
on E (without shadowing contributions from T g and <*) requires couples 
of points differing from each other only with respect to r # . It is 
inpossible to match with the data available such a constraint. Then 
points with almost the same value of Tg, but with different albedos 
were selected from the data-set. The data for these points are given 
in Table 7. 

Table 7. Experimental input data corresponding to days and sites with 
T s and o-data similar to those in table 6, but with high 
aerodynamic resistances, r . The last column again computed 
with eq. (18) 



h 

5* ® 
w 

T a 

(K) 

a 

r 

3 -1 

(day. cm ) 

E 

(mm. day ) 

E‘(T ,a) 

8 -1 

(mm. day ) 

1 

308.2 

308 

0.477 

2.33 

4.4 

2.09 

2 

307.2 

309.4 

0.422 

2.68 

2.3 

2.58 

3 

307.6 

309.8 

0.442 

2.43 

2.7 

2.4 


It must be emphasized that the term E (eq. 21) is independent 

~a 

from the formulation given by eq. (18). Therefore taking care of the 

sign, the variation <5 E must be added to the values of the column 

r a 

headed E'(T ,a). According to eq. (21) the variations 6 E between 
8 r a 

each separate point in Table 6 in combination with all the points of 

Table 7 were calculated. Results are shown in Table 8. From Table 8 

it can be concluded that the calculation of ^E requires only 

approximate values of r . When r can be assigned to about 0.1 day. cm"' 1 

a a _ | 

for a rough surface, and about 2 day. cm for a smooth one, 6 E is 

r a 

determined sufficiently accurate. In Table 8 the row-index applies to 
the rows of Table 6 and the column-index to the rows of Table 7. 
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Table 8. Values of 6 E (mm. day ) evaluated from Tables 6 and 7. The 

r a 

sign of 6 8 is for decreasing values of r 



1 

2 

3 

Average 

I 

1.65 

1.5 

1.5 

1.55 

2 

0.9 

0.9 

0.9 

0.9 

3 

1.0 

1.0 

1.0 

1.0 

4 

3.5 

3.3 

3.05 

3.3 


For each row averages can be taken and can be used for estimation 

of corrections for differences in r . 

a 

Similar calculations as performed for Table 5 yield the results 

shown in Table 9. Improvement in the estimation of E is evident when 

comparing the standard deviations of E and E-{e'(T )+<$ E|. It is 

w s,a r a J 

important to recall that, unless Table 4, only those data were used 
with D T.,a < °- 


Table 9. Comparison between experimental and calculated values. Eq. 
(18) corrected according to eq. (21) 


E-function 

Mean 

0 


(mm. day ) 

(mm. day ') 

E 

3.6 

1.25 

E'(T a)+6 E 
s » *a 

3.3 

1.4 

E-(E’(T s ,a)+6 r e] 

A 

0.3 

0.9 


d. Effectiveness of air temperature 

The results of Table 9 show that calculated values are rather 
close the measured data. However standard deviation is still too high. 
When looking at the data it seems that a systematic underestimation 
occurred during the period from 2 September to 8 September, 1978. This 
period was characterized by very high values of air temperature; the 
mean air temperature amounted to 35.8°C, while being 27.4°C over the 
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remaining data. Thus it was decided Co evaluate an additional correct-* 
ion from eq. (8), using 6T fl ■ 8,4°c. Accordingly d^E * 1.8 net. day 1 
was found. As for 6^ E the correction d^E is independent from the 
previous ones and must be added to E'(T ,a). Applying such a correct- 
ion yielded the results s..cwn in Table 10. 


Table 10. Comparison between experimental E and calculated (eq. 18) 

evaporation data with the corrections d E (eq. 21) and 

r a 

d E (eq. 8) being included 


E-function 

Mean 

(mm. day *) 

0 

(ran. day - * ) 

E 

3.6 

1.25 

E* (T ,a)+d_ E+d T E 
s r a 13 

3.7 

1.28 

E-lE’(T s ,a) + d r E + d Ta E] 

8 4 

-ti.l 

0.5 


e. Summary of the results 

Experimental values of actual evaporation rate, from different 

surfaces and periods, were compared with values calculated from eq. 

(18). Further corrections for the ae-odynamic resistance r as obtained 

a 

from eq. (21) and for the air temperature T a as obtained from eq. (8) 
were applied (Table 10). The accuracy of the fit is also shown in 
Fig. 5, where measured versus calculated values are plotted. An inter- 
val equal to 20 of E^jV (T g ,d) + E + d Ta is also shown in Fig, 5. 

In the evaluation of the corrections, the effect of atmospheric 
instability was not taken into account. However such an effect may be 
very iiiportant as far as period 2/9-8/9 is concerned. In those days 
surface temperatures even higher than 60°C were observed. In these 
conditions strong buoyancy phenomena develop. For these reasons the 
evaluation of r g by eq. (19) under unstable conditions might be unre- 
liable. 
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actual ♦voporation (cole ) 

Fig. 5. Measured versus calculated (eq. 18) evaporation rates. Correct- 
ions included relate to surface aerodynamic resistance and air 
temperature. A scatter equal to twice the standard deviation 
is shown (broken lines) 


IV. SUMMARY AND CONCLUSIONS 

A simple procedure to evaluate actual evaporation has been derived 
by a geometrical representation of the energy balance equation, ex- 
pressed as its first order Taylor's expansion. The role of the physical 
variables involved in the energy balance equation has been analysed. 

A linear relationship between actual evaporation and surface tempera- 
ture was calculated and compared with experimental data. A bi-linear 
relationship between actual evaporation, surface temperature and 
albedo was also determined. Comparison with experimental data yielded 
promising results. 

Corrections were applied as related to differences in aerodynamic 
resistance and air temperature for each day -experiment. The slope of 
the relationships between actual evaporation and surface characteristics 
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was shown Co be strictly tied v: aerodynamic resistance. 

Ho relevant effect of the indatermination of the soil heat flux 
in desert dry soils was found. 

After corrections were applied agreeaent between calculated and 
experimental evaporation data was shown to be good. 

The foraulas derived are oriented to applications involving remo- 
tely sensed surface teoperature and albedo. The required basic data 
are: aetaoro logical data, height and kind of soil coverage, aean value 
of the evaporation rata. Thee-* data allow the calculation of first- 
stage evaporation formulas. Short tern experiments over different ex- 
treme situations were shown to be more useful then long-lasting expe- 
riments over e few surfaces. 

■ * 

The linear functions derived can be used for a straight forward 
evaluation of actual evaporation from MSS and IRLS data. These data 
provide the required values of surface temperature end albedo for 
each single ares. The present procedure can be followed as a scheme, 
when evaluating evaporation losses from large inhomogeneous areas. 
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LIST OF USED SYMBOLS 


Symbol Interpretation 

a surface reflectivity 

E evaporation rate 


Units 


“ 2 a ' 

g.cm .day 


E(T g ,«) 

E f (T s ,a) 


e 

e' 

e* 



evaporation race calculated as an 
arbitrary function of T t a 

evaporation rate calculated as a 

linear function of T , a 

s 

surface emissivity 

air apparent emissivity 

ratio of molecular weight of water 
vapour to dry air 

soil heat flux 

height of roughness elements 
latent heat of vaporization of water 
atmospheric pressure 
thermal capacity of air 
surface aerodynamic resistance 
shortwave incoming radiation 
Stephan-Boltzmann constant 


air temperature 

T g surface temperature 

D E(x, t ... partial derivative of E with respect 
xi 1 % 

n' to x^. E is any function of 
(x j , ... x n ) 

6 E(x., finite variation of E with respect 

x i x ) 

n to x^. E is any function of 

(Xj, ... x n ) 


>2 

g.cm • 


diy 


I 


mm. day 


-1 


J.cm .day 


cm 



Pa 

J.cnf 3 .K -1 

day. cm * 

. "2 . 'I 

.T.em .day 

-2 - ; 
J.cm .day .K 

K 

K 
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